arXiv:1506.03678vl [physics.flu-dyn] llJun2015 


Under consideration for publication in J. Fluid Mech. 


1 


Diffusive dynamics and stochastic models of 
turbulent axisymmetric wakes 

G. Rigas f , A. S. Morgans, R. D. Brackston, J. F. Morrison 

Department of Aeronautics, Imperial College London, London SW7 2AZ, UK 

(Received ?; revised ?; accepted ?. - To be entered by editorial office) 

A modelling methodology to reproduce the experimental measurements of a turbulent 
flow under the presence of symmetry is presented. The flow is a three-dimensional wake 
generated by an axisymmetric body. We show that the dynamics of the turbulent wake- 
flow can be assimilated by a nonlinear two-dimensional Langevin equation, the determin¬ 
istic part of which accounts for the broken symmetries which occur at the laminar and 
transitional regimes at low Reynolds numbers and the stochastic part of which accounts 
for the turbulent fluctuations. Comparison between theoretical and experimental results 
allows the extraction of the model parameters. 
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1. Introduction 

Turbulent flows are ubiquitous in natural phenomena and engineering applications 
( Popd 2Q00h therefore a mathematically tractable description of them is desirable for 
their prediction and control. At low Reynolds numbers, corresponding to laminar regimes, 
bifurcation theory has aided understanding of the dynamic behaviour of fluid flows in 
the presence of symmetry: under certain assumptions, the infinite-dimensional nonlinear 
fluid system governed by the Navier-Stokes equations can be re duced to a simple and 
finite-dimensional system close to the threshold of bifurcation (jCrawford fe Knobloch 
ll99l[lHolme~ et aZ.lf2012h . However, departure from the laminar regimes and the critical 
bifurcating points renders the flow chaotic and finally turbulent, increasing the order of 
the system and the complexity for a mathematical description of it. 

Here we analyse the turbulent wake generated by an axisymmetric bluff body at high 
Reynolds numbers. This type of flow u ndergoes a steady bifu r cation followed by an un¬ 
steady one at low Reynolds numbers ( Bohorauez et al. 20 111 : iBurv fe Jardinl 2012h the 
dynamics of which can be accurately captured using simple weakly nonlinear models 
( Fabre et ali 20081 : Mcliga et al. 20091 ) prior to the emergence of chaos and turbulent 
behaviour. During the transitional regime, continuous spatial and temporal symmetries 
are spontaneously broken. Although bifurcations break the symmetries en route to tur¬ 
bulence, fully developed turbulence is known to restor e the possible symmetries in a 
statistical sense, at very high Reynolds numbers ( Frisch 19961) . 

In this paper we show that the dynamics of the turbulent wake can be approximated 
by a simple nonlinear Langevin equation, the deterministic part of which accounts for the 
broken symmetries which occur in the laminar and transitional regime, and the stochastic 
part for the turbulent fluctuations. 

The paper first presents an overview of axisymmetric wake-flows and models in the 
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laminar (S3 and turbulent (S3 regimes, before presenting the stochastic model in 
i}2j The model is validated using experimental data ()J3|), and the results are presented 
in fjT] 


1.1. Laminar regime and modelling 

The laminar and linearly stable axisymmetric wake loses spatial rotational symmetry, TZ^, 
in the azimuthal direction (rotation of angle 7r around any radial ax is passing through 


the centre of the body) due to a supercritical pitchfork bifu rcation (|Fabre et gZ.1 12008 : 


Meliga et qZ.ll2009t Bohorauez et aZj|201ll : [Bobinski et aZ.ll2Q14 1. the normal form of which 
reads 


x = ax 


Ax|x| 2 (1.1) 

(a, A £ R, A < 0). Above the threshold of instability (a > 0), Eq. (11.11) is associated 
with symmetry breaking since the steady state solutions are not invariant under the 
x —> —x symmetry. This model can be directly obtained from the Navier-Sto kes equations 
throu gh a weakly nonlinear expansion around the critical bifurcating point ( Meliga et all 
2009) and has been used extensively for the description of laminar flows undergoing 
supercritical bifurcations ( Drazin fe Rcidj 2004ll . In the context of the weakly nonlinear 
analysis, x = x(t) describes the complex amplitude evolution of the bifurcated state 
(global mode). 

Increasing the Reynolds number, a subsequent Hopf bifurcation is responsible for the 
loss of continuous temporal symmetry and the emergence of period ic shedding of vortice s 
before the turbulent axisymmetric wake reaches a chaotic regime ( Burv & Jardinl 2012 1. 


1.2. Turbulent regime 

At high Reynolds numbers, a large number of studies have shown that the bifurcating 
states observed in the laminar wakes persist and manifest as large-scale structures. For 
the axisymmetric wake, it was recently shown that the steady bifurcated state (respon¬ 
sible for the loss of rotational symmetry) and the unsteady bifurcated state (responsible 
for the vortex shedding) persist at high Reyn olds numbers and fully turbulent regimes 
( Rigas et al. 2014 : Grandemangc et al ] 20L4). Similar behaviour has been observed in 
wakes with ot her types of symmetries (see for example the rectilinear three-dimensional 
Ahmed body ( Grandemange et al. l l2012l.l201.llB . 

Here we extend the laminar weakly nonlinear modelling approach to the fully turbu¬ 
lent regime, and we show that the derived model captures the dynamic evolution of the 
turbulent three-dimensional wake. This is achieved by modelling the effect of the tur¬ 
bulent fluctuations acting on the deterministic dynamics of the system, reminiscent of 
the transitional instabilities observed in the laminar regimes, as stochastic forcing. The 
stochastic modelling approach of the turbulent dynamics has been successfully applied 
to various turbulent flows exhibitin g global and large-scale dynamics including the tur¬ 
bulent Rayleigh-Benard convection ( Sreeniva san et al. 2002: Brown fe Ahlersll2007l) and 
the turbulent von Karman swirling flow ( de La Torre fc Burguetef 2007 1. Hereafter we 
demonstrate our modelling approach considering the dynamics associated with the first 
steady bifurcation given by Eq. m- 


2. The stochastic model 

We consider the deterministic system given by Eq. m and apply a stochastic mod¬ 
elling approach for the turbulent background fluctuations. In Cartesian coordinates, 
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Figure 1. Stationary probability density function (left) and potential (right) of the nonlinear 

Langevin model. 



x = ( x,y ), the deterministic system with independent additive white noise becomes 

x = ax + Xx(x 2 + y 2 ) + cr^fft), 

y = ay + Xy(x 2 +y 2 ) + cr^y(t), (2.1) 


where £ accounts for the random forcing of turbulence and a 2 the variance of it. Here, 
a = cr x = cr y and the stochastic process in Eq. (EUD is rotationally symmetric. We 
transform the system from ( x , y) —1 (r, <)>), where r is the radial dista nce from the ce ntre 
(amplitude) and (j) the angle (phase), using the Ito interpretation (jGardine r ,1985 ). In 
polar variables, the Langevin system given by Eq. becomes 


r = ar + X r 3 + 


2 r 


er£ r 



( 2 . 2 ) 


where r = \Jx 2 + y 2 and cf> = tan _1 (y/x). Notice that in polar variables the radial 
component is independent of the angular position. 

The stationary probability density function (PDF) of the above system can be found 
from the steady state Fokker-Planck equation and is given by 


P s (r , cf) 


C exp 


U(r, <t>) 

K 


(2.3) 


where C is a normalisation constant, K = cr 2 /2 is the noise intensity (diffusivity) and U 
the potential. The potential U is 


U (r, <t>) 


2 

ar 

~y 


Ar 4 

-h K In r 

4 


(2.4) 


The stationary PDF for the angle is uniform, J 0 °° P s (r, f>)dr = ^. In the case of a rota¬ 
tionally symmetric experimental setup and inflow conditions, the drifts and diffusivities 
are independent of <p. The stationary PDF, Fig. [T] peaks around the minimum of the 
Mexican-hat shaped potential. 


3. Experimental set-up 

The turbulent wake is generated by an axisymmetric bluff-body with a blunt trailing 
edge, a schematic of which is shown in Fig. [2] The body diameter D , is 196.5 mm and 
the length-to-diameter ra tio L/D , is 6.48. M ore details for the expe rimental setup and 
procedure can be found in R.igas et al. ( 2014 ) and lOxlade et al .1 (1 2015 ). Experiments were 
performed at a constant free stream velocity, t/oo = 15 m/s. The Reynolds numbers based 
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Figure 2. Schematic representation of the model, (a) Side view; (b) rear view. Pressure mea¬ 
surements are taken on the base of the body from 64 static pressure tappings, based on which 
the Centre-of-Pressure location is calculated. 

on diameter and boundary layer momentum thickness at separation are Rep = 1-88 x 10 5 
and Reg = 2050 respectively. The boundary layer separation and the ensuing wake are 
turbulent. 

Pressure measurements are taken on the base of the body from 64 static taps equally 
spaced on a uniform polar grid with Sr = 0.056.D and S(f> = 45° in the radial and 
azimuthal directions, respectively. The sampling frequency was 225 Hz. A total of 19,200 
s of data was acquired over sixteen independent experiments, providing approximately 
2000 independent measurements with a 95% uncertainty of approximately 0.45% and 1% 
in time average and rms pressure respectively. 

4. Results 

Let x of Eq. ed represent the Centre-of-Pressure (CoP) coordinates. The CoP is 
used as a macroscopic variable quantifying the symmetry of the turbulent wake. Here 
it is calculated from the space-averaged pressure, defined on the Cartesian coordinate 
system of the base, x' = (x',y'), and non-dimensionalised with the body diameter, as 

x(t) =Mm //«>*'"• 

where A is the area of the base of the body. A zero CoP value corresponds to an 7Z n - 
symmetric flow, w hereas departure from this value to an increased asymmetry of the flow 
( Rigas et ~al\ 12014ft - 

The temporal evolution and PDF of the radial and azimuthal location of the CoP 
for the highly turbulent regime are shown in Fig. [3] Highly erratic motion with time 
is observed in both components. Statistically, the wake spends most of the time in a 
non-zero radial location, indicating broken 1Z^ symmetry. Due to the uniform PDF in 
the azimuthal location, rotational symmetry is recovered in the long time average. For 
the laminar regime and after the first bifurcation, which is described by Eq. ED, the 
stable equilibrium corresponds to a non-zero CoP, the angle of which is determined from 
the initial conditions and is unique. During the turbulent regime, the laminar stable 
equilibrium point explores a continuum of states around its mean radial value and uni¬ 
formly distributed in the angular direction: hence, the wake explores an infinite number 
of metastable states restoring the lost lZ n symmetry. 

Insight into the random dynamics of the turbulent wake is provided from the calcu¬ 
lation of the time-averaged mean-square displacement (MSD), defined as ([Ax(r)] 2 ) = 
((x(f + t) — x(t)) 2 }. The MSD of the angular and radial components is plotted in Fig. [4] 
for different sampling times. In the azimuthal direction, the MSD increases linearly with 
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Figure 3. A time series and probability density of the Centre-of-Pressure: angular (upper) and 
radial (lower) position. Symbols in the probability density (right) correspond to experimental 
data and solid lines to theoretical model predictions. Dashed line: mean radial value indicating 
the broken 7symmetry. 


a A K = a 2 12 

3.81 -5604 0.0028 

Table 1. Coefficient values of the Langevin model obtained experimentally at ReD = 1.88xl0 5 . 
The units of the coefficients are s _1 since lengths are non-dimensionalised with the body diam¬ 
eter. 


time, ([A</>(t)] 2 ) oc r, consistent with free diffusive motion (lEinsteinl 119051) . In the ra¬ 
dial direction, the linear relation holds only for short time scales below a threshold 
t)] 2 ) oc t, t < t s , and reaches a saturation plateau at larger time scales, 
lim T _ ) . 00 ([Ar(r)] 2 ) = H 2 . 

The above results for the CoP are consistent with the predictions of the model given 
by Eq. ([2T2]) . The coefficients of the Langevin equation a, A and K , are obtained from 
the experimental MSD and PDF and are given in Tabic [T] The slope of the radial MSD 
relation is directly correlated with the diffusion coefficient K, which is obtained through 
linear fitting, ([Ar(r)] 2 ) = 2 Kt, t < t s . Knowing the diffusion coefficient K, the model 
coefficients a, A are uniquely defined from Eq. m and here are obtained through least- 
square fitting. 

In Fig.[3l we compare the measured probability density function with the one predicted 
by the model from Eq. (12.31) . In Fig. 0] we compare the MSD of the experimental results 
with the numerically calculated MSD from the model. Direct numerical integration of 
Eq. (12.21) was performed using an Eulcr-Maruyama scheme. We notice that the dynamics 
of the CoP can be described over all the time scales from the Langevin model. 

A qualitative explanation of the above results is provided by the potential well shown 
in Fig. [0 The turbulent wake, the state of which is quantified by the CoP, meanders in 
the Mexican-hat-shaped potential and explores an infinite number of states through a 
random walk (diffusive motion). Specifically, in the azimuthal direction, which determines 
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T (S) 

Figure 4. Mean square displacement as a function of the time interval r. (a) Angular and (b) 
radial components of the CoP location. Symbols: experimental data. Solid line: model. Thick 
dashed lines: power-law fits. Insets show power spectral density of the experimental CoP com¬ 
ponents. An exponent of -2, consistent with free diffusive motion, is obtained for the azimuthal 
and radial components. At low frequencies, the PSD of the radial components is constant due 
to the spatial confinement imposed by the potential well, in accordance with the MSD results. 


the orientation of the wake, it can explore freely any azimuthal location resulting in 
unbounded reorientations. In the radial direction, the motion is restricted due to spatial 
constraints (confinement imposed by the potential well) resulting in a constant MSD at 
large timescales. 

The diffusive dynamics of the turbulent wake are also depicted in the power spectral 
density of the CoP, plotted as inset in Fig. [4j The spectral density is closely related 
to the MSD: in general, for a power law behaviour of the MSD, ([Ax(r)] 2 ) oc r K , the 
asymptotic form of the power spectral density is $(/) oc /~( 1+K ). In accordance with 
the linear increase of the MSD with time (k = 1), an exponent of -2 is observed in 
the spectrum of the angular component consistent with Brownian motion. A similar 
decay is observed for the radial component when / > l/t s . However, at low frequencies 
corresponding to / —► 0, or equivalently large timescales, it reaches a plateau and levels 
off, in accordance with the MSD measurements for r —> oo. 
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Figure 5. Reorientation characteristics: Joint PDF of angular variation and radial position of 
the CoP as measured from the experiment. Inset shows normalised conditional probabilities: 
data (symbols) collapse on a Gaussian distribution (line), as described by the model. 


The dynamics associated with the radial and angular motion of the CoP have been 
analysed independently so far. Now, we analyse the coupled dynamics predicted by the 
model. Specifically the coupling arises as an inverse relationship between <j) = Acj>/At 
and r for the angular component, as described by Eq. (12.211 . The model suggests that the 
conditional PDF of <j> for a given r follows a Gaussian distribution with variance inversely 
proportional to r 2 , that is 




T\ff2/K 




(4.1) 


Indications of the inverse relationship can be found in the time series of the CoP in 
Fig. [3l One observes that for small values of the radial component, abrupt changes of 
the azimuthal orientation occur. The inverse relation is further validated from the joint 
PDF of the angular variation and radial location, p(r, A<fr), shown in Fig. 0 where large 
reorientations are more probable at small radii. The conditional PDFs of the reorienta¬ 
tions A (f>, p(Acf>\r) = p(r, Atf>)/P(r) collapse to a zero-mean Gaussian distribution with 
variance er 2 , when scaled by 1 jr and plotted against r<j>. In the limit of small radii val¬ 
ues, r —> 0, angular rotations with an almost uniform PDF are observed in the joint 
PDF. This situation of small r corresponds to large variance and therefore a very broad 
Gaussian that, over the finite measurable windo w of — 7r to 7r, will appear uniform. This 
is analogous to the cessation events observed bv lBrown et al. ( 20051 ) during the reorien¬ 
tation of the large-scales of Rayleigh-Benard convection and can be understood without 
specifying arbitrary thresholds for a small r. 


5. Conclusions 

The physical picture that can be drawn for the turbulent axisymmetric wake based 
on the above is as follows. The laminar large scale structures, associated with spatially 
broken symmetries, persist at high Reynolds numbers. In the turbulent regime, these 
structures undergo diffusive motion (random walk) in a two-dimensional Mexican-hat- 
shaped potential well (see Fig. [l} restoring statistically the broken symmetries. In this 
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paper we have shown that this behaviour can be captured by a simple nonlinear model 
consisting of two coupled stochastic differential equations. The deterministic part of 
the model accounts for the broken symmetries observed in the laminar regime, and the 
stochastic part models in a phenomenological sense the turbulent fluctuations acting on 
the large scale structures. The unsteady bifurcation, responsible for the vortex shedding 
(limit cycle), has not been considered in the present analysis. This can be done by adding 
an ex tra equation accounting for the temp oral broken symmetries due to a Hopf bifurca¬ 
tion ( Fabre et «/J 120081 : Meliga et al. 2009h : experimentally, the quasi-periodic nature of 
the vortex shedding appears in the PSD of the CoP as peak at ~ 16 Hz corresponding 
to a Strouhal number of ~ 0.2 based on D and U ac . 

The diffusive dynamics of the large scale structures presented here and in iRigas et~al\ 
( 20141) show close similarities with the diffusive dyna mics of the large-scale circulation ob¬ 


serve d in the turbulent Rayleigh-Benard convecti on (jBrown et al 


2005: Brown fc Ahlers 


2006) and the turbulent vo n Karman swirling flow ( de La Torre fe Burguetel2007t Berhanu et al 
2009 1” 


Petrelis et al. 2009f ). However, this similarity is not surprising if the symmetries 


of the problem are accounted for: spatial rotational symmetry which breaks en route to 
turbulence. 

In summary, we have shown that turbulent dynamics can be described by determinis¬ 
tic equations derived from symmetry arguments with stochastic forcing terms that give 
rise to turbulent dynamics. The modelling process we have presented here should be ap¬ 
plicable to other turbulent flows, provided that their specific symmetries are taken into 
account. 


This work was supported by the Engineering and Physical Sciences Research Council 
(grant number EP/I005684). 
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